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ABSTRACT 

We present new calculations of the evolving UV background out to the epoch 
of cosmological reionization and make predictions for the amount of GeV gamma- 
ray attenuation by electron-positron pair production. Our results are based on re- 
cent semi-analytic models of galaxy formation, which provide predictions of the dust- 
extinguished UV radiation field due to starlight, and empirical estimates of the con- 
tribution due to quasars. We account for the reprocessing of ionizing photons by the 
intergalactic medium. We test whether our models can reproduce estimates of the ion- 
izing background at high redshift from flux decrement analysis and proximity effect 
measurements from quasar spectra, and identify a range of models that can satisfy 
these constraints. Pair-production against soft diffuse photons leads to a spectral cut- 
off feature for gamma rays observed between 10 and 100 GeV. This cutoff varies with 
redshift and the assumed star formation and quasar evolution models. We find only 
negligible amounts of absorption for gamma rays observed below 10 GeV for any emis- 
sion redshift. With observations of high-redshift sources in sufficient numbers by the 
Fermi Gamma-ray Space Telescope and new ground-based instruments it should be 
possible to constrain the extragalactic background light in the UV and optical portion 
of the spectrum. 

Key words: intergalactic medium - gamma rays: bursts - gamma rays: theory - 
cosmology: theory - diffuse radiation - ultraviolet: galaxies 



1 INTRODUCTION 

Interactions between photons via electron-positron pair pro- 
duction can have a substantial effect on the observed spectra 
of extragalactic gamma-ray sources. This process, which can 
occur when the required threshold energy of twice the elec- 
tron mass is present in the centre-of-m ass frame, removes 
gamm a rays en route to the observer (|Gould fc Schredeil 
1 19671 ). and provides a link between high-energy astrophysics 
and the extragalactic background light (EBL) - the inte- 
grated luminosity of the universe at UV, optical, and IR 
frequencies. 

Attempts to measure the present-day EBL through 
absolute photometry in the optical and IR portion of 
the spectrum are hampered by the bright foregrounds 
of the Milky Way and zodiacal light from dust in 



E-mail:rgilmore@physics. ucsc.edu 



the solar system, as well as calibration uncertainties 
dHauser fc Dwekl 1200 ll ). Measurements in the near- and 
far-IR are available from th e DIRBE and FIRAS expert 
ment on the COBE sat ellite (IWright fc Reesell2000l: IWrighi 



Cambrcsv ct al. 2001; Lcvenson, Wright & Johnso^ 
2007: Hauser et al...l99a: Lagache et al.„ 2000. ; Wright , 200i : 



Fixsen et al.lll99M). A tentative determination of the optic al 



EBL was made inlBernstein. Free dman fc Madord (|2002al lbh 
with revision in Bernstein (2007). Number counts of ob- 
served galaxies provide a robust lower limit to the EBL, but 
the degree to which these measurements converge can be 
controversial. The analysis of the integr ated galaxy counts 
in seven optical and near-IR bands by iMadau fc Pozz ettil 
(2000) showed that flux from direct starlight converged 
to a fairly low level, below that claimed in the sev- 
eral DIRBE detections of near-IR flux. Recent experi- 
ments such as the Galaxy Evolution Explorer (GALEX) 
satellite, ISO, and Spitzer have provided counts data in 
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a variety of non-optical wavelengths, placing lower lim- 
its on the local EBL at most wavelengths of interest. 
A more thorou gh discussion of current EBL constr aints 
can be found in Gilmore. Primack fc Somervillel (|2009l ) and 
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A number of techniques have been employed to build 
cosmological models of the emission of light from galax- 
ies at UV to far-IR wavelengths. These include interpola- 
tion and extrapolat ion of cosmological observables such as 
star formation rate (iKneiske. Mannheim fc Hartmannll2002l : 
iKneiske et al.ll2004l : iRazzaaue. Dermer fc Fin kell2009l). and 
lumin osity functions (Franccschini, Rodighie ro fc Vaccaril 
I2OO8I ). Another method evolves backward in time the lo- 
cal galaxy population, usually by assuming that the lu- 
minosity density changes with redshift as a pow er law at 
all wavelengths (|Stecker. Malkan fc ScuUvl 120061 ). A third 
class of techniques models the galaxy population forward 
in time, beginning with cosmological initial condition s 
l|Primack et al.ll200ll : IPrimack. Bullock fc Somervilldlioosh . 

Our goal in this paper is to build a suite of models 
of the evolving background light produced by stars and 
quasars, with a focus on the optical-UV background out to 
high redshift. A preli minary report on this work appeared in 
iGilmore et al.l l|2008l ). We calculate the reprocessing of ion- 
izing radiation by the intergal actic medium (IGM) usin g the 
radiative transfer code CUBA (|Haardt fc Madau|[l99^ . and 
use the observed ionization state of the IGM to constrain 
our models. The contribution of starlight to the EBL is 
predicted by recent semi-analyti c model s (SAM s) of galaxy 
formation, de scribed in S onicrv ille et al.l (|2008|). These have 
been used bv I Gilmore et al.l (see also IPrimack et al] 

I2OO8I ) to make predictions for the evolving EBL. Most previ- 
ous modeling attempts, including our recent work with this 
new SAM, have focused on observational data in the optical 
and IR. These are the wavelengths most relevant to obser- 
vations of relatively nearby (2 < 0.5) blazars with ground- 
based instruments, which until recently have typically fea- 
tured energy thresholds above ~ 150 GeV. With the recent 
launch of Fermi with its Large Area Telescope (LAT) sen- 
sitivity range of 20 MeV to 300 GeV, as well as the advent 
of new ground-based experiments such as MAGIC-II with 
energy thresholds < 100 GeV, it is now important to make 
theoretical predictions of the UV background at ionizing and 
non- ionizing wavelengths out to high redshift. This paper is 
an attempt to specifically target absorption in this region of 
the gamma-ray spectrum. 

Understanding the absorption that occurs for gamma 
rays observed between 1 and 100 GeV is an uncertain un- 
dertaking due to the lack of sensitive observations of the 
EBL at the corresponding UV wavelengths. Moreover, the 
dechning opacities for gamma rays in this region means that 
sources are likely to be visible out to large redshift. Evolu- 
tion of the background must be taken into account when 
calculating absorption for all but the nearest blazars, and 
at high redshifts the EBL can have a spectral energy dis- 
tribution (SED) much different than observed locally. The 
most distant object with confirmed redshift that has been 
detected at VHE e nergies is currently the flat-spectrum ra- 
dio quasar 3C279 (|Albert et al.1 12OO8I ) a,t z = 0.536. This 
object was observed at energies between 90 and 500 GeV, 
with a steep spectrum that was likely due in part to EBL 
absorption. In the 10 to 100 GeV energy decade that is now 



being probed by Fermi and upgraded ground experiments, 
the characteristic redshift at which the EBL becomes op- 
tically thick to pair-production is expected to increase to 
redshifts of several. 

A small number of calculations have been performed 
that speciflcally addressed the question of the gamma-ray 
absorption by the UV background. In iMadau fc Phinnevl 
(1996), two different models of star-formation, based on dif- 
ferent assumptions about the B-band normalization, were 
used to predict gamma-ray opacities from 10 to 200 GeV, 
with propagation of ionizing photons through the IGM taken 
into account. This work suggested that the universe becomes 
optically thick at a few tens of GeV for gamma rays emitted 
at 2 ~ 2. A secon d work which focused on the UV back- 
ground, [o3 (|200ll ). argued that the absorption by ionizing 
photons was negligible, and that < 20 GeV observed gamma 
rays would only be significantly attenuated at higher red- 
shifts, where they would interact with photons below the 
Lyman limit. Lya photons were found to be a significant 
component of the UV flux. This paper also explored the pos- 
sibility of using Fermi to detect an evolving blazar attenu- 
ation edge, which would probe hi gh-redshift star-formation . 
Finally, the background model of lSalamon fc Steckerl (Il998l ) 
targeted absorption of 10-500 GeV gamma rays, and used 
an estimate of high-redshift star-formation based on evolu- 
tion seen in damped Lya systems. This work also included 
a UV contribution from quasars. 

At energies below the Lyman limit, lower bounds on 
galaxy emissivity exist from number counts by GALEX 
(Xu_et al. 2005) and Huhhle S pace Telescope (HST), as well 
as ba lloon-based experiments ([Gardner. Brown fc FergusonI 
I2OOOI ). Such experiments are subject to systematic errors 
in completeness and photometric measurement of apparent 
magnitude, and can only test the background out to mod- 
erate redshift. At higher redshifts we no longer have mea- 
surements that directly connect to the EBL, such as direct 
number counts and absolute photometry, and uncertainties 
and possible biases in cosmological measurements such as lu- 
minosity functions and star-formation rate density become 
increasingly problematic. 

Measurements of the ionization state of the IGM can 
provide constraints on ionizing flux. At redshifts higher than 
the 'breakthrough redshift' k. 1.6, the universe is optically 
thick to Lyman continuum photons, and ionizing fields be- 
come local, with a mean free path that decrea ses rapidly at 
larger redshifts (|Madau. Haardt"fc" Reedl 19991 ) ■ while below 
this redshift the mean free path becomes longer than the 
horizon length. Studies of the opacity of Lya and other red- 
shifted absorption lines place constraints on the emission of 
UV photons by probing the neutral fraction, and therefore 
the balance between photoionizations and recombinations 
(|Haehnelt et al.ll200ll : iMadau et al]|2004al ). As these lines 
are affected by the local radiation, they provide information 
about sou rces existing at approxi mately the redshift of the 
absorber (|Haardt fc Madau|[l99^ ). 

Two methods of determining the ionization state 
of the IGM include the proximity effect, in which one 
searches for the decrease in Lya emiss io ns near an AGN 
JPairAglio. Wisotzki fc Worseckl I2OO8I : iLiske fc Willigeil 
l200ll) . and flux decrement analysis, which utilizes hy- 
drodynamic simulations to model the distribution of 
Lya absorption along the line of sight to an AGN (e.g. 
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iBolton et all |2005| ). The line-of-sight proximity effect 
utilizes the decrease in absorption lines in the vicinity of a 
quasar, compared to farther away along the line of sight, due 
to increased ionization fraction. As the quasar has a known 
UV luminosity, the deficit of absorption in this region can be 
used to estimate the background; a larger change indicates 
a lower background flux. As quasars do not reside in typical 
cosmological environments, a number of potential biases 
exist. Quasars tend to be found in overdense environments, 
which can lead to overesti mates of the background flux by 
as much as a factor of 3 (|Loeb fc EisensteinI Il995l ) . Time 
variation in luminosity on the time-scale of photoionization, 
typically ~ 10* years, will also tend to bias results towards 
a high background, as quasars tend to be selected in their 
brigh test phases l|Schirber. Miralda-Escude fc McDonald! 
I2OO4I '). It is also now recognized that using broad emission 
lines such as Lya tends to lead to underesti mated redshifts 
and therefore higher quasar luminosity l| Richards et al.l 
I2OO2I ). This may have been a problem in many determi- 
nations of the proximity effect. The assumed cosmological 
model also affects the resulting background inferred by 
these measurements. The second met hod mentioned, th e 
less-direct flux decrement technique l|Rauch et al.l Il997l ). 
is not without its own potential biases; it relies on correct 
cosmological parameters and knowledge of the quasar's 
unabsorbed continuum level, a problem at high redshift 
where absorption is strong. Newer attempts to correct 
for the biases in proxi mity effects measurements, such as 
iDall'Aglio et al.| l|2008h , have found lower values for the 
ionizing background flux that are more consistent with the 
flux decrement technique. 

Observations of the Lya forest can also provide clues 
about the types of sources producing the ionizing back- 
ground, which in our model include star-forming galax- 
ies and quasars. The quasar luminosity function (LF) has 
been measured by large-scale surve ys such as the Two- 
degree Field (2dF) (iBovle et al.ll2000l) and the Sloan Digital 
Sky Survey (SPSS ) l|jiang et al.ll2006l : [Richards et aLll2006l : 
ICroom et al.|[2003 l. and data are also available at a variety 
of frequ encies from experiments such as XMM, Chandra, and 
Spitz er (iBarger et al.ll2005l : iMatute et al]|2006l : iBrown et all 
I2QO6I ). The hydrogen of the intergalactic mediu m (IGM) is 
know n to be fully ionized below a redshift of ~6 l|Fan et al.l 
l2006h . Photons above the Lyman limit are responsible for 
reionizing the universe and maintaining it in a highly ion- 
ized state. The relative contributions of star-forming galax- 
ies and AGN to this process are not fully understood, but 
there is evidence that quasars are a sub-dominant com- 
ponent at this epoch. The decline of the quasar luminos- 
ity function observed beyond redshift three constrains the 
quasar contribut i on to the ionizing b ackground to be ~ 10~^ 
l|Fan et al.ll200ll : iMadau et ai]|l999l ). unless there is an un- 
expected steep upturn in the quasar luminosity function at 
low lu minosities. A new approach bv ISrbinovskv fc Wvithj 
l|2007l ) utilizing semi-analytic modeling sets limits on the 
quasar contribution to ionizing radiation of 1.4 to 14.5 
per cent at z = 5.7, and studies of the soft X-ray back- 
round also constrain this fraction to be subdominant 
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Increased quasar emission is believed to be responsible 
for He II reionization, which as tracked by He II Lya ab- 
sorption takes place at a lower redshift than hydrogen, z ^ 



3 (e.g. iBolton et al.ll2005l ). The shape of the ionizing back- 
ground therefore evolves in redshift, with a hardening of the 
spectrum that is indicative of an increased contribution from 
quasars. The degree to which AGN dominate the UV back- 
ground at the time of He reionization is a debated issue, 
with some suggestions that stars and AGN pro vide roughly 
equa l contributions to th e background at z ~ 3 l|Kriss et al.l 
I2OOII : ISmette et al.l l2002V One of the major sources of un- 
certainty in this transition lies in the unresolvable faint en d 
of the AGN luminosity function (ISchirber fc Bullockllioosi '). 
The ratio between hydrogen and helium ionization fractions, 
particularly H I (13.6 eV) and He II (54.4 eV), can be used to 
measure the slope of the total UV spectrum in this regime. 
The decrease in the optical depth of He II indicates that the 
harder radia tion from quasars increases with t ime between 
z = 5 and 3 (|Shull et al.ll2004l : iFan eral]|2006l) . 

Aside from quasars, the known dominant sources of 
UV radiation are short-lived massive stars, mainly of O- 
and B-type, which closely trace the star formation rate den- 
sity. Estimating the ionizing contribution from star-forming 
galaxies directly is complicated by the fact that only a small 
fraction fesc of this radiation escapes from galaxies due to 
neutral gas and dust in the interstellar medium, as we will 
discuss in Section 2.4. 

Because of the uncertain nature and evolution of sources 
of ionizing radiation, in this paper we consider four mod- 
els that attempt to span a realistic range of assumptions. 
In Section 2 we discuss the inputs to our model, including 
a short description of our semi-analytic models, assumed 
quasar luminosity density, and radiative transfer code. In 
Section 3, we introduce our four UV background models and 
present results, including the evolving background radiation 
and comparisons with Lya forest measurements. The main 
results of the paper, gamma-ray opacities, are presented in 
Section 4, with a discussion following in Section 5. Our four 
models are summarized in Table [T] and their successes and 
failures in accounting for the data are summarized in Table 

m 



2 MODELING 

To calculate the evolving UV background we have used pre- 
dictions of the UV luminosity density from galaxies, as pro- 
vided by recent semi-analytic models (SAMs) of galaxy for- 
mation, combined with estimates of quasar emissivity. The 
combined emissivities from galaxies and quasars are inte- 
grated over redshift to find the evolving background flux. 
Photons from these sources at energies above the Lyman 
limit can be absorbed and reradiated by the IGM; we cal- 
culate the effect of these processes using the CUBA radiative 
transfer code. In this section we discuss the semi-analytic 
models and radiative transfer code used in calculating the 
background, and show some general results. 



2.1 The semi-analytic models 

The semi-analytic models used in calculating 
the EBL in t h is wo rk are described in detail in 
ISomerville et aL 1 (|2008l, SOS), and are based on the 
models described in | Somerville fc Primackl (|l999D and 
ISomerville. Primack fc Faberl l|200lh . with several new 
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updates and capabilities; we provide a very brief summary 
in the followins paragrap hs. I nterested readers sho uld also 
see lPrimack et al.l (2008) and iGilmore et al.l (j2009l ). which 
focused on the optical and IR EBL resulting from galaxies 
in this model. 

The backbone of the semi-analytic models are dark mat- 
ter "merger trees", which describe the cosmological assem- 
bly history of dark matter haloes as they build up over time 
through mergers and accretion of diffuse material. These 
merger trees may either be constructed via Monte Carlo 
techniques based on the Extended Press-Schechter theory 
(as in this work), or extracted from N-body simulations. 
We compute the rate at which gas can cool via atomic 
processes and be accreted onto the central galaxy within 
the potential well of each of these collapsed and virialized 
DM haloes. Cold gas is then converted into stars in both 
a "quiescent" and starburst modes. The "quiescent" mode 
takes place in isolated discs and is modeled according to 
an e mpirical prescription based on the Schmidt-Kennicutt 
law (jKennicutt et al.1ll998l ) . The starburst mode is triggered 
during mergers, with gas being rapidly converted into stars 
on time-scales determined from hydrodynamic simulations 
of galaxy mergers. Feedback from supernovae and massive 
stars can reheat cold gas and expel it from the galaxy. Chem- 
ical evolution is modeled using a simple instantaneous recy- 
cling approximation with the effective yield treated as a free 
parameter. 

The most recent models also treat the growth of su- 
permassive black holes within galactic nuclei and the im- 
pact of the released energy on galaxies and their environ- 
ment. Seed black holes, perhaps originating from the re- 
mants of very massive Pop III stars, are planted in the top- 
level haloes in the merger trees, and can grow by two ac- 
cretion mechanisms. Galaxy mergers trigger rapid (Edding- 
ton limited) "bright mode" accretion and correspond obser- 
vationally to classical X-ray or optically luminous quasars 
and AGN. Bright mode AGN activity drives winds that can 
expel cold gas from the galaxy. Accretion of hot gas from 
the hot halo feeds low-level (Bondi) "radio mode" accre- 
tion, which is associated with the production of giant radio 
jets. These jets are assumed to be capable of heating the 
surrounding hot gas halo, offsetting or even quenching cool- 
ing flows. Although the radio mode feedback in particular 
has been shown to play a key role in reproducing important 
galaxy properties in these models, such as galaxy colour bi- 
modality and luminosity functions, the predictions for ob- 
servable properties of quasars and AGN in these models has 
not yet been thoroughly tested. Therefore, in this work, we 
instead add in the contribution to the background radiation 
due to quasars empirically. 

The star formation and chemical enrichment histories 
for each galaxy are convolved with stellar population models 
to produce syntheti c spectral energy distribu tions (SEDs). 
We make use of the lBruzual fc CharlotI (|2003l ) models with 
a Chabrier stellar initial mass function. 

Our models include a simple estimate of the attenuation 
of starlight due to dust, based on a two-component model 
similar to that proposed bv , Char lot fc Fall (2000,). One com- 
ponent is the diffuse "cirrus" dust in the disc and another is 
associated with the dense 'birth clouds' surrounding young 
star forming regions. The l/-band, face-on extinction optical 
depth of the diffuse dust is given by 



Tdust.O ^cold n^cold 

" ' 

where Tdust,o is a free parameter, Zcoid is the metallicity of 
the cold gas, mcoid is the mass of the cold gas in the disc, 
and rgas is the radius of the cold gas disc. To compute the 
actual extinction we assign each galaxy a random inclina- 
tion and use a standard 'slab' model. Additionally, stars 
younger than lO'^ yr are enshrouded in a cloud of dust with 
optical depth tbc,v = Mbctv.o, where pec = 3. Finally, 
to extend the extinction correction to other wavebands, we 
assume a Galactic attenuation curve (|Cardelli et al.|[l989l ) 
for the diffuse dust component and a power-law extinction 
curve A\ oc (A/5500A)", with n = 0.7, for the birth clouds. 
The free parameters are adjusted to reproduce the observed 
ratios of far-UV to far-IR luminosity in nearby galaxies. 

We consider two different choices of parameters for our 
semi-analytic models, which differ primarily in the assumed 
cosmological parameters. The free parameters that control 
galaxy formation in each model are then tuned to match lo- 
cal galaxy observations, as described in SOS (the actual val- 
ues of the parameters for both models are also given in SOS, 
Table 2). The 'fiducial' model adopts a concordance cosmol- 
ogy with D,m = 0.3, Qa = 0.7, h = 0.70, and as = 0.90. Our 
'low' model adopts the best fit values from WMAP3 for these 
parameters, with i},n = 0.23S3, Qa ~ 0.7617, h = 0.732, and 
as ~ 0.761. The most relevant difference in this work is the 
value of the power spectrum normalization as- The lower 
normalization of the primordial power spectrum in the 'low' 
model leads to delayed structure formation and decreased 
luminosity densities at high redshifts (see SOS). Adopting 
the best-fitting parameters from the more recent analysis o f 
WMAP5, which favored as = 0.S2 (|Komatsu et al.H2009l ). 
gives nearly identical results to our WMAP3 (low) models. 

The predictions for attenuated and unattenuated emis- 
sivity from our fiducial semi-analytic model are shown for 
several redshifts in Fig. \T\ Predictions from the 'low' model 
are qualitatively similar, although because of the delayed 
star formation, galaxies tend to have higher gas surface den- 
sities, and therefore higher dust opacities and larger atten- 
uation values are predicted. 

In our companion E BL paper l|Gilmore et al.ll2009l . see 
also lPrimack et al]|2008l ). we compare the predictions of our 
SAMs with a broad range of data, including local luminosity 
functions, optical and IR luminosity density, and number 
counts in a variety of bands from the UV to far-IR. In a 
planned future paper we will make a careful comparison of 
these models with the galaxy population at high redshift. 
Here we show just a few representative results demonstrating 
that our models are doing reasonably well at reproducing 
the UV properties of local and distant galaxies. In Fig.[2]we 
compare our model predictions with galaxy number counts 
in two UV bands, using data from the GALEX satellite and 
other experiments. This provides a test of the low-redshift 
normalization of our model in the UV range. GALEX has 
surveyed the sky t hese bands and provided data down to 
magnitude ~ 23.5 l|Xu et al.lboosi '). At fainter magnitudes, 
there are measurements from the STIS instrument on the 
HST ((Gardner et a l. 2000), albeit with large uncertainty due 
to poor statistics. Populations of brighter objects have also 
been probed by the FOCA balloon-borne UV telescope, and 
counts from this instrument at 2000 A have typically yielded 
higher numbers than GALEX after wavelength correction. 
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Figure 1. The emissivity due to galaxies predicted by our fiducial 
galaxy-formation model at a number of redshifts, normalized to 1 
Mq/jt. Dotted curves show the emission predicted in the absence 
of dust extinction. Vertical lines indicate the ionization energies 
of H I and He II at 912 and 228 A. 



possibly due to differences in calibration. Our models show 
good agreement with the data at 2310 A, but are a bit higher 
than the GALEX observations at 1530 A, though they are 
not in disagreement with the FOCA data. 

Recent data from a variety of instruments has con- 
strained the UV luminosity density out to high redshift. 
In Fig. |3] we have compared the UV emissivity in galaxies 
from our models against data at a rest frame emission wave- 
length of approximately 1500 A. We find that data from the 
GALEX-VVDS, GOODS, and deep HST ACS imaging aU 
agree reasonably well with the UV evolution of our fiducial 
model. 



2.2 Star formation 

The star-formation rate density (SFRD) as a function of red- 
shift in the 'fiducial' and 'low' models are shown in Fig. U 
compared with observational estimates of star formation 
density at various redshifts, all of which has been converted 
to a Chabrier initial mass function (IMF). At z < 1, both 
of our models are in good agr eement with the observational 
compilation of lHopkin^ (|2004| ) , while at 1 < z < 2 they tend 
to skirt the lower envelope of observational values. However, 
there are still large discrepancies between SFR estimates 
from different indicators and different data sets at these red- 
shifts, in part due to the increase d fraction of sta r formation 
in heavily obscured systems (e.g. lHopkins|[2007l ). where the 
correction for dust obscuration is uncertain. At z > 2, the 
SFRD in the 'low' model declines fairly steeply, while in 
the fiducial model the SFRD remains nearly constant from 
2 < 2 < 5 and then declines more gradually. As discussed 
above and in SOS, this is because of the lower normalization 




15 20 25 30 

AB magnitude 




15 20 25 30 

AB magnitude 



Figure 2. Number counts in the GALEX 1530 A (upper) and 
2310 A bands (lower). The solid black line shows the fiducial 
model, and dashed blue shows the low model. Note that the low 
model has counts equal to or slightly greater than the fiducial 
model at some magnitudes; this is due to differing amounts of 
dust extinction at low redshift b etween the two models. Data 
are from GALEX jXu et al.ll2005l . green squares), STIS on HST 
jGardner et al.l l2000l. purple asterisks), and the balloon-borne 
FOC A experiment lllglesias- P aramo et al.l 12004 iMilliard et al.l 
19921. red stars and open pentagons respectively). Following 
Xu et al. lEooi), all counts have been converted to the GALEX 
bands by assuming a UV spectral slope of -0.8. 
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Figure 3. The emissivity at 1500 A as a function of redshift 
in our models. As previously, the solid black line is the fidu- 
cial model, and the dashed blue line shows the lo w model. The 
blue c ircle at redshift 0.1 is GALEX data from IWvder et al.l 
||2005| ') and the purple stars are measur ements using GALEX 
and other data llSc hminovich et al]|2005h . The red squares are 
from GOODS llDahlen et al.ll2007h . the red circles are determi- 
nations from ground-bas ed observations meddv et al 1 l2008h . and 
the green stars are from lBouwens et al.l | |2007^ . 



of the power spectrum and reduced small scale power in this 
model, which delays the formation of structure. 

Above redshift four, observational estimates of global 
star-formation rates diverge, and different measurements 
can disagree by as much as an order of magnitude. Studies 
of UV luminosity functions of dropout galaxies by Bouwens 
and collaborators ( Bouwens et al. 2008, 2 007 ) find relatively 
low values for the global SFR, with a monotonic decrease 
above redshift four. Higher rates have been found by other 
authors, including those who have derived s tar-formation 
histo ry from detections of gamma-ray bursts l|Yiiksel et al.l 
l2008h . These studies suggest a much higher rate of star for- 
mation which does not decrease significantly until z > 6. 
This may be due in part to the fact that the Bouwens et al. 
data points that we report here were obtained by integrat- 
ing the UV luminosity function down to a luminosity corre- 
sponding to 0.04 times the observed value of L« at redshift 
three. Other authors make different choices for the lower 
limit of integration, and the relatively steep slope of the UV 
LF at these redshifts implies that this can make a significant 
difference. The SAM predictions shown include the star for- 
mation in all galaxies (down to the mass resolution of the 
simulation) . 

Our 'low' semi-analytic model, based on the lower de- 
termination of as in WMAP3, produces a star-formation 
rate that is lower than most of the data points at mid and 
high redshift, and reproduces the rapid fall-off in star forma- 
tion indicated by the Bouwens points. Our 'fiducial' model, 
based on WMAPl, does a better job of matching the higher 



star-formation rates seen in other dropout analyses, as well 
as the data from gamma-ray bursts. 

Although not shown here, another way to constrain the 
star formation history of the Universe is via the build-up 
of stellar mass in the form of long-lived stars. As shown in 
SOS, the fiducial model predicts a stellar mass density that 
is higher than observational estimates by a factor of two 
at 2 = 1 and a factor of three at z = 2, while the 'low' 
(WMAP3) model produces good agreement with the stellar 
mass assembly history. As discussed in SOS, this discrepancy 
between the observed stellar mass assembly history and the 
star formation hist ory has been pointed out in a number of 
recent papers (e.g. lDave|[200^ : iFardal et al.ll2007D . and one 
possible explanation is that the stellar IMF was more top- 
heavy in the past (so that more UV photons were produced 
per total unit mass of star formation). However, it is also still 
possible that the discrepancy is simply due to observational 
uncertainties in the estimates of star formation rates and/or 
stellar masses. 

In addition to the star-formation histories predicted 
self-consistently in our fiducial and low semi-analytic mod- 
els, we consider an additional ad hoc high-peaked fiducial 
form for the SF history above z = 3. This is not a semi- 
analytic model; it is simply a functional form that was cho- 
sen to be consistent with the highest observational deter- 
minations of the star- formation rate. We then utilize the 
same redshift-dependent dust extinction factors as the fidu- 
cial model. We include this case to illustrate the predictions 
for gamma-ray attenuation for an extreme model with the 
maximum plausible UV background at high redshift. How- 
ever, we note that as the fiducial model already produces an 
integrated stellar mass density in excess of that observed at 
high redshift, the high-peaked model is strongly disfavoured 
by these observations. 



2.3 Radiative transfer 

Ionizing photons from galaxies and quasars which escape 
into the intergalactic medium (IGM) are processed by neu- 
tral hydrogen and neutral and singly-ionized helium which 
resides in Lya forest clouds (LAC) and thicker Lyman- 
limit systems (LLS), defined here as having column densities 
> lO^'^'^ cm~^ . This has a strong effect on the spectrum and 
intensity of the average background field. The propagation 
of ionizing flux through the IGM in our models is calculated 
using an updated version of t he CUBA code. An earlier ver- 
sion of CUBA is described in iHaar dt fc Madau' (2001') and 
is based on the theory of iHaardt fc Madau (1996). Here we 
briefly summarize some of the main ideas and formalism 
from these papers. 

The effect of residual neutral gas on the ionizing radia- 
tion field can be described in general terms by the radiative 
transfer equation: 



, a ad. 



-3 — J — ckJ + - — e. 

a 47r 



(2) 



Here Jif) is the intensity of the radiation field for frequency 
f, e{v) is the emissivity, a is the cosmological scale factor, 
c is the speed of fight, and k is the continuum absorption 
coefficient. This equation accounts for both the redshifting 
of photons to lower energies, and absorption by neutral gas. 
Quasars and star-forming galaxies contribute to e{p) in our 
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Figure 4. The global star formation rate density predicted by our 
models, compared with a compilation of observational data. The 
solid black and dashed blue curves show the SFRD history of our 
fiducial and low models, respectively. The black dash-dot curve 
which diverges from the fiducial curve above redshift three is the 
'high-peaked' form which we discuss in the text. The red s quare s 
at lower redshift are from the compilati on of iHopkinj l|2004h 
The purple stars are from observations by iBouwens et al.l teOOsl . 
l2007t ) of dropout-selected galaxies. For these we show the dust- 
corrected results from integrating the luminosity functions down 
to a value of 0.04 L* at 2 = 3; it is possible that fainter objects 
provide an additional contribution. The magenta squares at red- 
shift 4.5 and 6 show inferred star forma tion rates from gamma-ray 
burst observations l lYiiksel et al.|[2"008t) . The green c ircle is based 
on ob servations of Lyman-break galaxies at z ~ 5 llVerma et al.l 
|2007|). and the orange triangle is an upper limit from VLT data 
llMannucci et al.l2007l). The blue s quares are results from the Sub- 
aru Deep Field llOuchi et "alll2004l ') . All data have been corrected 
for extinction (by the authors) and converted to a Chabrier IMF. 



model, along with the diffuse reemission of absorption sys- 
tems. Lyman absorbers are taken to have a distribution that 
can be described in terms of power laws in column density 
and redshift 



(X n: 



^(1 



(3) 



dNnidz 
with parameters 

7 = 0.16 (LAC, < 2 < L4) 
7 = 3.0 (LAC, 1.4 < z) 
7 = 1.5 (LLS, all redshifts) 

used in these models. A distribution with slope -1.5 in col- 
umn density has been shown to des cribe absorbers over 
a wide range in Nhi (|Hu et al.l Il995l ). and the slopes for 
redshi ft evolution are based on observational determina- 
tions (iKim et al ] | 19971 : IStengl er-Larrea et al ] | 19951 : iBechtoidI 
I1994I ). The effective optical depth from the Lya forest ab- 
sorption in this distribution is shown in Fig. [5] Note that 
our LyQ optical depth does not follow the upturn seen at 



Figure 5. Effective optical depth as a function of redshift from 
our assumed absorption cloud distributi on. Data from quasar 
spectra are shown at (5 < z < 6.5) from lFan et al.l (2006); here 
the blue points are averaged Gunn-Peterson measurements, and 
the red, cyan, and green symbols are Lya, /3, and 7 measurements 
of the hi ghest-redshift i ndivid ual objects. Values at lower redshifts 
Schave et al.l l(20o3) (green circles) and lDall'Aglio et al.l 



are from 



ll2008h (red stars). 



z ~ 6, where a rapid rise in absorption may signal that our 
assumption of a uniform UV background is no longer valid. 

It should be noted that the exact form of the column dis- 
tribution function is not critical, as it is the integrated value 
of this parameter from which we derive the effective optical 
depth and therefore the average background. The effective 
depth is dominated by clouds with opacity near unity. Us- 
ing a power-law form simplifies the integration process, and 
speeds up our computation with little loss in accuracy. We 
do caution readers that the choice of distribution function 
can have a large effect on He II absorption, and in turn on 
the background above 54.4 eV which we report in Fig. 1101 
The background at these energies is not expected to affect 
our gamma-ray attenuation significantly because the photon 
density at these high energies is so low. 

Lyman systems reradiate a fraction of the absorbed 
light via recombination radiation. Our code accounts for the 
contribution of H I recombinations to the UV flux via free- 
bound, Lya, and two-photon continuum emission. For the 
latter two, only the non-ionizing background is affected. The 
total proper volume emissivity from IGM clouds from radi- 
ation released for a particular mode can be quantified as 



e{u,z) = hvft {v) Wabs{z) '~{z,v) ^ 

atot al 



(4) 



where Ui is the fraction of recombinations leading to the 
particular mode, which has probability fiii') of creating a 
photon of energy v. In the case of Lyman-a emission this is 
simply a delta function at the line energy, and fo r the contin- 
uum distributions descriptions can be found in lOsterbrockl 
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l| 19891 ). The remaining functions contain the details of emis- 
sion and absorption from absorption systems 

4TrJ{u', z'' 



hv' 



dNuidz 



-Wab. 



:{v')dv' 



Pem{ly, NHl)dNE 



(5) 



(6) 



The first of these quantities is the rate of ionizations by 
the background field J(i^', z). Here Wabsiy') encodes the in- 
formation about the optical depth for photons of a given 
energy, and thus the probability of being absorbed. Equa- 
tion [6] is the integral over absorption systems, multiplied by 
Pem.{i', Nhi), which is the probability of a photon of given 
energy escaping from a cloud after emission. The code does 
not include the contribution from s awtooth mod ulation due 
to H I and He H Lyman resonances (|Madau fc H aardt 200^). 

We do not include radiation from recombinations to 
He I, as neutral helium is small in number density com- 
pared to both H I and He H (though that may not be the 
case at redshifts o n the verge of reionization) . As shown in 
iHaardt fc Madaul (1996), thermal coUisional effects can pro- 
vide a sizable fraction (20 to 30 per cent) of the H I emission 
by Lya and two-photon processes; these are not accounted 
for in our code. Collisions between He H atoms are never 
significant, as there is insufficient thermal energy to excite 
these modes. 



2.4 Ionizing escape fraction from galaxies 

The fraction of ionizing photons produced that escape from 
galaxies into the IGM is a free parameter in our model. This 
parameter is poorly constrained, with observations and sim- 
ulations giving widely different and sometimes conflicting 
results. In the literature, the escape fraction may be defined 
in a couple ways. The absolute escape fraction is simply the 
fraction of radiation at wavelengths just shortward of 912 
A which escapes the dust and neutral hydrogen in a galaxy. 
This definition is most relevant for the purposes of modeling 
the ionizing background. What is actually measured in ob- 
servations is the relative escape fraction, where the ionizing 
flux is compared to a non-ionizing wavelength, often 1500 
A. As described below, we use a relative definition which 
separates the amounts of attenuation due to dust and H I 
(Equation [7]) . Ionizing radiation from quasars is not atten- 
uated in escaping the host galaxy in our model. 

Direct detection of escaping UV radiation has only been 
successful in a handful of individual cases. As two of these 
detections have bee n for galaxies at Z'^ S with large escape 
fractions measured l|Shaplev et al the fact that many 

low redshift attempts to find ionizing radiation have failed 
with low upper bounds may suggest evolution in this quan- 
tity between redshifts one and three. Rather firm upper lim- 
its on escape fraction from direct detection effo rts exist for 
lower redshift galax i es (se e the compilation in ISiana et al.l 
120071 ). ISteidel et all |2to3) reported ionizing flux from 29 
stacked galaxies at z ~ 3.4, at a level indicating little or no 
attenuation. However, this result suffered from a selection 
bias, as the Lyman-break galaxies used were chosen from 
the bluest quartile of the population. 

While observations have mainly determined upper lim- 
its on the ionizing escape fractions, some authors have used 



the ionization state of the IGM to derive lower limits. Ion- 
ization rates inferred from Lya forest data and reasonable 
extrapolations of source number to faint luminosity can re- 
quire a high escape fraction. Values of ^ 20 per cent above 
redshi ft 5 were found to be needed in iBolton fc Haehneld 
(|2007l ). ISrbinovskv fc Wvithd i 20081 ) have found that con- 
straints on the escape fraction from 5.5 < z < 6.0 from 
N-body simulations require a global minimum of 5 per cent 
to match Lya data, with a higher fraction needed in the 
event that star formation in galaxies in smaller haloes is 
suppressed . 

Recently, detailed Adaptive Mesh Refinement (AMR) 
N-body hydrodyna mical simulations of high-re dshift galax- 
ies (3 < z < 9) by iGnedin. Kravtsov fc ChenI ((2008 ) have 
found low escape fractions of 1-3 per cent, without strong 
evolution in redshift or dependence on galaxy properties. 
This work found that most escaping ionizing radiation orig- 
inated from stars in a thin shell at the outside of the H I disc. 
Smaller galaxies have less escaping radiation due to the fact 
that their H I discs are thicker relative to the distribution of 
young stars. Dust is found to have little effect on the escape 
of Lyman-continuum radiation, as the unobscured minority 
of stars that provide most of the escaping ionizing radia- 
tion have essentially no attenuation due to dust, while stars 
that have translucent (r ~ 1) optical depths due to dust are 
generally completely obscured by HI. Another analysis has 
been undertake n using a smoothed particle hydrodynamic s 
(SPH) code by iRazoumov fc Sommer-Larseni (|2007l . I2006D 
and has found evolving escape fractions, with feac = 6 to 10 
per cent at z = 3.6 decreasing to 1 to 2 per cent at z — 2.4. 
This simulation did not include the effects of dust. 

Our semi-analytic models predict the emissivity from 
star-forming galaxies down to a minimum rest-frame wave- 
length of 100 A. Lyman-continuum photons are attenuated 
by a factor /esc, hi; this determines the absorption of pho- 
tons shortward of 912 A and is a non-evolving input to our 
radiative transfer code. This parameter is defined as the fol- 
lowing ratio of intrinsic and observable luminosities at 912 
and 1500 A 



(-£'gi2/il50o) e 



n 



/esc/i50C 



(7) 



where /esc is the absolute attenuation factor from both dust 
and H I for ionizing photons near the Lyman limit, and 
/i500 is the factor from dust alone at 1500 A. As dust ab- 
sorption is an evolving effect included in the semi-analytic 
model, total absorption /ese for ionizing photons escaping 
from galaxies must be interpreted as the product of /cse,Hi 
and the average /isoo for a particular redshift in the model. 
The average value of /isoo in our model is higher at high 
redshift, and is higher in the low model than in the fiducial 
model at high redshift. Typical values of /i^qq, as seen in 
the difference between attenuated and unattenuated curves 
in Fig. [ll vary from about 6.8 at 2 = 2 to 11.2 at z — 5 
for our fi ducial model, a nd 9. 7 to 14.2 in the low model. As 
noted in iGnedin et al.l (|200^) , taking the total absorption 
to be a product of the dust and H I factors may not be 
physically realistic, as the two components are may not be 
distributed in the same way within the galaxy, however it is 
sufficient for the purpose of calculating emissivities, as we 
are dealing only with global quantities. 
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2.5 Quasar emissivity 

Quasar input to our model is accomplished using an as- 
sumed UV luminosity density, which determines the out- 
put at all energies via a fixed spectral form. We have used 
the quasar luminosity fu nctions at 912 A determined by 
iHopkins. Richards fc HernauistI (|2007l . HRH07) which are 
based on a large observational data set and spectral and ob- 
scurational modeling. This work found that, with the appro- 
priate corrections for obscuration, a single bolometric func- 
tion could match data in each band. Both the bright- and 
faint-end slopes of the LF are argued to become more shal- 
low at higher redshift, indicating a increasing contribution 
from bright, unobscured AGN, which dominate the total 
AGN luminosity above redshift ~ 2. 

Another approach to modeling t he quasar contribu- 
tion to the background was presented in ISchirber fc BuUockl 
l|2003l . SB03). This work estimated the evolution of the 
unobserved faint-end of the quasar luminosity function at 
high redshift using observational constraints on the ionizing 
background. They found very different results depending on 
whether they used constraints on the ionizing background 
from the quasar proximity efi^ect or from the Lyman- a flux 
decrement method. Assuming the broken power law univer- 
sal quasar spectrum we present below, the HRH07 luminos- 
ity density evolution at 912 A is similar to Model 'A' of 
SB03, which produced the lowest fluxes and was found to 
be consistent with the flux decrement data when combined 
with a substantial contribution from star-forming galaxies. 
The highest derived flux, arising in Model 'C', was suffi- 
cient to produce ionizing photons at the level suggested by 
proximity effect measurements. While we recognize that it is 
probable that many of these proximity effect measurements 
overestimate the background due to aforementioned biases, 
we will use the SB03 Model 'C' for the purposes of creating 
and analyzing an extremely quasar-dominated background 
model. 

Luminosity densities at the H I and He II Lyman limits 
for each of these models are shown in Fig. |6l We have renor- 
malized the original model proposed by SB03 by a factor of 
0.8 to better match the HRH07 results below redshift 2.3. 

Quasars are assumed to have a spectrum which can 
be modeled as a broken power law in F^, = dF/du. In the 
ex treme-UV, w e have adopted the hard spectrum suggested 
bv lTelfer etaP (|2002l l. from observations with the Faint Ob- 
ject Spectrograph on HST of quasars at a wide variety of 
non-local redshifts {z > 0.33). The values we assume are as 
follows 



F, oc 

with indices 

Q = 0.4 
a = 1.3 
a = 0.2 
Q = 0.5 
Q = 1.57 
a = 0.9 



(8) 



(12 ^im < A) 
(1 /^m < A < 12 /im) 
(500 nm < A < 1 /im) 
(120 nm < A < 500 nm) 

(A < 120 nm) 
(Soft: X-rays, > 500 eV) 



The contribution of quasars to the background is only non- 
negligible in the extreme-UV (A < 120nm). At other wave- 
lengths we use slopes to roughly fit typical x-ray spectra of 
Seyfert 1 galaxies, and the optical and IR emission of quasars 
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Figure 6. Quasar lumi nosity density for t he iHopkins et al.l 
(2007) (solid black) and ISchirber fc Bullock feOOSl) Model 'C 
(broken green) models at 912 (upper lines) and 228 A (lower 
lines). The latter has been multiplied by a factor of 0.8 to better 
match the observations at low redshift. 



seen in SPSS and other observations l|Vanden Berk et al.l 
I2OOII : ISanders et a"Llll989l ) . 



3 COSMOLOGICAL MODELS AND 
RESULTING UV BACKGROUND 

To calculate the evolving background radiation field, we 
have combined three estimates of the star-formation rate 
density with two possible forms for the quasar luminosity 
de nsity. The two quasa r models are the 'realistic' estimate 
of lHopk ins et al] l|2007h from a large body of observational 
data, and a higher 'extreme' model motivated by proximity 
effect measurements for 3 < 2; < 5 from SB03. Our four 
models are summarized in Table [T] 



3.1 Lya forest constraints 

Inferred ionization rates and column density measurements 
from the Lya forests of quasar spectra provide us with an 
independent measurement of the UV background through 
its ionizing properties, which we can compare with the IGM 
state computed by our radiative transfer code. In Fig. [T] we 
compare the ionization rate (in terms of r_i2, the average 
rate per hydrogen ion with units of 10~^^ s"'^) with data 
from both the quasar proximity effect and the flux decre- 
ment in Lya forest measurements. As discussed in the In- 
troduction, these two techniques have tended to give dis- 
parate values for r_i2. For the flducial and low models with 
the HRH07 QSO LF, we show ionization rates with several 
values for the escape fraction of ionizing radiation from the 
galactic H I disc. With a moderate H I escape fraction of 0.1 
to 0.2, our flducial model is able to reproduce the level of 
ionizing background detected by most determinations using 
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Table 1. The background models considered in this work. The 
second and third columns show the star-formation histories and 
quasar luminosity densities used as inputs in each model. The 
escape fraction in the last column refers to the values used in cal- 
culating the background flux and optical depth to gamma rays in 
Figs. llOl and im respectively. Our star formation history scenarios 
are discussed at the beg inning of this section . 'HRH07' refers to 
the best-fitting m odel of iHopkin s et aP 1I2OO7I '). and 'SB03 Model 
C to the model in lSchirber &: B ullock (2 0031) . We have multiplied 
the latter by a factor of 0.8 to better match the observed quasar 
luminosity density at low redshift. The escape fraction refers to 
the attenuation of ionizing photons from star-forming galaxies by 
neutral hydrogen; attenuation by dust is included intrinsically in 
our semi-analytic model. 



Model SFR Density Quasar Luminosity /esc, HI 

1 Fiducial HRH07 0.1 

2 Low HRH07 0.2 

3 Fid. High-peaked HRH07 0.1 

4 Fiducial SB03 Model C 0.02 



1 I — I — I — I — I — I — I — I — I— pi — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — p 




flux decrement techniques. Including attenuation by dust, 
this corresponds to a total ionizing escape fraction of ~ 1 to 
3 per cent, consistent with upper limits f rom observations as 
well as values suggested by simulat ions (jCnedin et al ] |2008l : 
iRazoumov fc Sommer-Larsenll2007l ). With the low model, a 
higher H I escape fraction of ~ 0.5 is necessary to match the 
highest redshift points (« > 5), due to the rapidly declining 
star- formation rate at high redshift. This escape fraction is 
higher than suggested by some authors, but there are no 
direct constraints on escape fractions at such high redshifts. 

Based on the quality of fits for our different models to 
these flux decrement data, we have chosen escape fractions 
of 0.1 and 0.2 for our fiducial and low models, respectively, 
as reasonable values to use in calculating the background 
and pair-production opacity. Both of these models predict 
ionization rates which decline above redshift 2.5; this is due 
both to the shape of the star-formation history and the in- 
creasing opacity of the IGM with redshift. 

The flux decrement calculati o ns fro m the l argest quasar 
sampl es, iFaucher-Giguere et al.l l|2008l) and iBolton et al.l 
lj2005l ). find ionization rates that are essentially flat from 
z = 2 out to z = 4, albeit with differing normalizations. Our 
model predictions are reasonably consistent with these ob- 
servational estimates, considering the uncertainties involved. 
The high-peaked model better reproduces the flatness of the 
ionization rate from z ~ 3 to 4, but still predicts too steep a 
rise with time from z ~ 2 to 3. The SF history predicted by 
our fiducial model could be made perfectly consistent with 
the Faucher-Gigu'ere et al. (2008) data by assuming an es- 
cape fraction that evolves from ~ 0.2 at z ~ 4 to 0.02 at 

2-2. 

The final scenario we e xamine uses the higher qu asar 
emissivities of Model 'C' in ISchirber k. Bullockl (120031 ). As 
the ionizing contribution from star-formation is subdomi- 
nant at all intermediate redshifts in this case, we have as- 
sumed a low escape fraction of 0.02. The ionizing flux in this 
model is capable of reproducing the highest measurements 
from the proximity effect. We have already mentioned sev- 



Figure 7. Ionization rate per hydrogen atom (with units of 
10~^^ s~^) in our four scenarios compared with data at a 
range of redshifts. Black solid lines: fiducial model with H I es- 
cape fractions from star-forming galaxies of 0.02, 0.1, and 0.2 
(bottom to top). Dashed blue lines: low star-formation model, 
with escape fractions 0.1, 0.2, and 0.5. Orange dash-dotted line: 
high-peaked star formation rate with escape fraction 0.1. These 
aforementioned models all use the quasar emissivity of HRH07. 
The long-dashed g reen line shows the fiducia l SFR model with 
quasar model C of ISchirber fc BuUockl | |2003|) . and escape frac- 
tion 0.02. Data points are divided into those obtained from 
fiux-decrement analysis (black) and those obtained via prox- 
imi ty effect near qua s ars (green ) . Ref e rences for the former 
arelBolton et al.l ll2005h:lFan et al] ll2006l):lFaucher-Giguere et al.l 



..20081): iKirkman et al] 1I2OO5I): iTytler et al.l ||2004|). and the 
latter i nclude IScott et al.l 1I2OOOI): ICooke. Espey fc Carswell 



poke. H/Spi 

1 19971) : iGialloneo. Fontana fc Madaul lll997^: ICristiani et al 
I I995I): IWilliger et al.l lll994l): ILu. Wolfe fc Turnshekl l|l99ll) 



I B a i tlik .Duncanfc'Ostrikerl 19881) . Some points have been 



shifted slightly for readability. 



eral known biases which may have artificially elevated these 
values, and this model should be considered an extreme pos- 
sibility. 

Another Lya forest measurement which can provide in- 
sight into the UV background is the relative abundance of 
H I and He II present in the IGM. This may be presented 
in terms of relative column densities A^(HeII)/A^(HI), or 
analogously as inverse ionization rates for these components 
FHi/rHoii; this is often referred to as the UV softness param- 
eter. In Fig. |8]we show how softness evolves with redshift 
for each of our background models. Our low, fiducial, and 
high-peaked star-formation densities with the HRH07 QSO 
LF are able to provide a reasonable match to observations 
when a moderate escape fraction is assumed. High escape 
fractions ^ 0.5, which are required for the 'low' SFR model 
to match ionization rates at high redshift, tend to overpre- 
dict softness. Our quasar-dominated model (SB03 model C) 
does not reproduce the trend of increasing softness in the 
background field with redshift, another factor which disfa- 
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Figure 8. The ratio of He II to H I column densities plotted 
against redshift. Higher values indicate a softer ionizing back- 
ground, with comparatively more ionizing photons available per 
hydrogen atom. Line types in this plot are the same as in Fig. [7] 
Data are fr om observations of H e II Lya systems bv lZheng et al.l 
1I2OO4I ') and lFechner et al.l 1 I2OO6I) . 



vors such a dominant contribution from faint quasars. The 
column density ratio in this case is not found to be sensitive 
to the H I escape fraction. 



3.2 The background flux 

The key result of this work is a prediction of the evolving UV 
background out to redshift 2; ~ 9, which has been calculated 
from our models for the total (stellar + quasar) emissivity 
e{i', z) combined with a calculation of the absorption and 
re-emission by IGM radiative transfer processes. Photons at 
non-ionizing wavelengths evolve passively, and are not at- 
tenuated significantly by any source. For these photons, the 
flux at a redshift z p and frequ ency uq in proper coordinates 
can be written as (|Peebles| [l99ai 



An 



dl (1 + zof 
dz (1 + ^)3 



e(i',z)e ^"f^dz, 



(9) 



where f = 1/0(1 + z) / (1 + zq) and dl/dz is the cosmological 
line element, defined as 

dl c 1 ^^^^ 



dz (1 + z)Hq ^!^„(l + z)3 + nA 

for a flat ACDM universe. Here Te// is the optical depth for 
ionizing radiation d ue to Poisson-distri buted HI absorption 
systems, defined as (|Madau et al.|[l999l ) 



re//(z/o, zo, z) 



dz 



dNn 



' dNuidz 



7(1 



(11) 



where r is the Lyman continuum optical depth through a 
given cloud, and N / {dNuidz) is the absorber distribution 
that we specified for our model in Eq. [3] For non-ionizing 



Figure 9. The predicted z=0 EBL spectrum from our semi- 
analytic models of gala xy formation from the UV to far-IR 
(see also iPrimack et al.||2008( ). We show results for the fiducial 
(black) and low (dashed blue) models, compared with experi- 
ments at a number of wave l ength s. Ot her lines show the E BL 
models of lFranceschini et al. 1 1|2008|) and lKneiske et al.l (120041 ') for 
comp arison. The blue-viole t triangles are results from STIS on 
HST llGardner et al.ll2000l). while the open magenta triangles 
are from GALEX l IXu et al] |2005|) . The green and red tria ngles 
are from the Hubble Deep Field l lMadau fc PozzettH |2000^ and 
Ultra Deep Field (Dolch, in preparation) respectively; the for- 
mer also includes near-IR grou nd based-data. O pen red trian- 
gles are from IRAC on Spitzer llFazio et al.ll2004) . The blue di- 
amonds are iBernsteinI ||2007| ) and other symbols in the near-IR 
are from several analyses of DIRBE data at 1.25 , 2.2, and 3.5 fim 
llLevenson fc Wright 200^; iLevenson et al.ll2007l: ICambresv et al.l 



I2OOII. iGoriian. Wright fc Charvl I2OO0I: IWright fc Reesell2000l). In 
the mid-IR, counts data is shown from ISO IIE lbaz et. a lll2002n 
and Spitzer jPapovich et al.ll2004l : [Charv et al.l2 004: Frave r et aU 
l2006l/lDole et al.11200 6'). In the far-IR, direct detection data is 
shown from DIRBE IWright..200 j. blue stars). 



radiation, the effective opacity is zero. At wavelengths well 
above the Lyman limit, the background at a redshift zo is 
determined by the total history of emission at higher red- 
shifts 2; > Zo- At ionizing wavelengths, the mean free path of 
photons is shorter than cosmological distances at redshifts 
greater than the breakthrough redshift, and therefore the 
ionizing background at zo > ztr is determined by the emis- 
sivity of approximately contemporary sources. 

In Fig. [9] we present the z = EBL fiux for our semi- 
analytic models, alongside two other recent EBL determina- 
tions using alternate methods. At low redshift, the effect of 
our high quasar and high-peaked star-formation models on 
the total background are negligible, and are not shown. We 
have shown the EBL calculated out to the far-IR in this plot, 
though results at wavelengths longer than the optical-near- 
IR peak at ~ l^m are not relevant to gamma-ray opacities 
in the high-redshift regime we discuss here. Dust reemis- 
sion of light in this model is based upon the IR templates 



12 R. C. Gilmore et al. 



100 




Figure 10. The history of the background flux, shown at the present day and three other redshifts. Intensities at the nonzero redshifts 
have been multiplied by the indicated factors (10"'^, 10~*, and 10~^ at z = 1.6, 3.2, 5.8, respectively) to separate the lines. Black solid 
line: fiducial model with H I escape fraction of 0.1. Dashed blue line: low star-formation model, with escape fraction 0.2. Orange dash- 
dotted line: high-peaked star formation rate with escape fraction 0.1. Green long-dashed: fiducial model with SB03 quasar contribution 
and escape fraction 0.02. At low redshift, only the first two models are shown, as the other models do not produce discernibly different 
levels of background at these times. We have also shown observational measurements of the background flux at 2 = in the UV, optical, 
and near-IR from Fig. [O] 



of iDevriendt fc Guiderdonil l|2000l ). Much more detailed IR 
background modeling, tak ing into account the latest Spitzer 
data, will be presented in iGilmore^ et al.l (|2009l ). 



models), and IGM reprocessing. The feature at ~300 A 
due to He II Lya. 



Understanding the evolution of the background field is 
critical to our calculation of the absorption of gamma rays, 
as the flux at high redshift can be both significantly higher 
than and of a different spectral shape than the local back- 
ground. In Fig.[TO]we show the background fiux for our four 
models at several redshifts, including z = where we have 
also shown a compilation of observable data, including es- 
timates from both absolute photometry and discrete source 
number counts. In each case, we assume the 'best-fitting' 
escape fractions of Table [T] which we have chosen based on 
the comparisons with Lya forest measurements. At all red- 
shifts, the background shows a sharp drop at the Lyman 
edge; this is a combined consequence of absorption in stel- 
lar atmospheres, H I in galaxies (quantified as /oscHI in our 



4 GAMMA-RAY ATTENUATION 

Gamma rays can interact with background photons if suf- 
ficient energy exists in the centre of mass to create an 
electron-positron pair with total mass 2me: 



(12) 



where Ei and E2 are the photon energies and 6 is the angle 
of incidence. A gamma-ray of energy can therefore pair- 
produce with background photons above a threshold energy 
of 



Et. 



E^{1 - cose) 



(13) 
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The cross-section for this process is (|Madau fc PhinnevI 
1 19961 ) 



16 



2/3(/3^ - 2) + (3 - r ) In 



1+P 
1-P 



where 



1 - 



2mic* 



£;i£2(l -cos6i)' 



(14) 



(15) 



and (Tt is the Thompson scattering cross section. The cross 
section is maximized for centre of mass energies of approx- 
imately twice the threshold energy 2meC^, and falls as ap- 
proximately as inverse energy for E ^ Eth- The likelihood 
of absorption is maximized for photons at about 4 times 
the absolute threshold energy, with one factor of 2 from a 
and another in going from 9 = tv ('head-on' configuration) 
to the most probable angle of interaction 9 « n/2. Gamma 
rays above 1 TeV are most attenuated by the near- and 
mid-IR range of the EBL, while those in the 300 GeV to 1 
TeV regime are sensitive to light in the near-IR and optical. 
Below 200 GeV it is mainly UV photons that have sufB- 
cient energy to cause the pair-production interaction. Below 
19 GeV only background photons with energies above the 
Lyman limit have sufficient energy to interact at any angle 
in the rest frame, and there is little attenuation. Note that 
these numbers refer to rest-frame energy of the gamma ray, 
which can be substantially higher than its observed energy 
due to redshifting. If the differential number density of back- 
ground photons at energy Ebg is n{Ebg, z), then the optical 
depth of attenuation for a photon of observed 

1 di 

t(Ej,zo) — - dz — / d(cos9) (1 — cos^) 

n oo 

X dEbgn{Ebg,z) a{E^{l + z),Ebg,9). (16) 

Where we have 

Emin = Eth + z) = 



E^{1 + z){l~ COS 9) 

and dl/dz is the cosmological line element defined in Equa- 
tion [TUl 

For each of our four models, we have calculated the op- 
tical depth of gamma rays at all relevant energies and red- 
shifts. As in our calculation of the background ffux above, 
we assume the H I escape fractions listed in Table [T] It 
should be emphasized that the choice of escape fraction has 
little effect on absorption of gamma rays at energies > 10 
GeV. We find, as argued in lOhl (|200lh . that the background 
field at energies above 1 Ry is negligible as a barrier to cos- 
mological gamma rays, and that significant optical depth 
above this energy is due to photons longward of the Lyman 
limit, where photon density increases dramatically in all of 
our models. While gamma rays are limited to interactions 
with background photons of an absolute minimum energy 
Eth = m1c^/E-y (with cos 6 = —1), redshifting places these 
gamma rays at higher energies at earlier epochs, where they 
can pair produce on the non-ionizing background. The in- 
crease of star-formation rate density by roughly an order of 
magnitude between present-day and peak rates means that 



gamma rays from high-redshift sources will tend to be at- 
tenuated most strongly at these early redshifts, where they 
have energies (1 -I- z) times higher than ai z — 0. 

In Fig. 111! the optical depth vs gamma-ray energy 
is shown for each model at various redshifts. These high- 
redshift results should be considered complementary to the 
our other calculations of EBL with these semi-analytic mod- 
els, which emphasized the absorption of > 100 GeV gamma 
rays at lower redshift. The effect of the UV background is to 
produce a relatively sharp and featureless cutoff in energy. 
At energies above 100 GeV, the effect of the EBL has often 
been quantified as a change in the spectral index of observed 
blazar spectra (e.g. lAharonian et al.|[200^ . due to the rela- 
tively fiat number density of EBL photons in the near and 
mid-IR. At lower energies, this approximation is not valid 
over any significant range in energy, due to the steepness 
of the cutoff that results in rapidly increasing numbers of 
photons with increasing wavelength in the UV. 

Our high-peaked SFRD and quasar-dominated models 
give absorption features that are similar, despite being very 
different in terms of the spectral form of the background ffux. 
While the emission from quasars produces a much higher 
ionizing background, the spectral cutoff at all redshifts we 
have explored is dominated by the photons longward of the 
Lyman limit. 

In Fig. 1121 we show the redshifts at which the universe 
becomes optically thick (r > 1) to gamma rays of a given 
energy for each of our models. From this plot, we can see in 
a general sense how background attenuation affects different 
energy regimes at different redshifts. The low model shows 
little change at redshifts higher than about 3, due to the 
rapid decline in star formation after this point. The high- 
peaked model has the most impact at high redshift, and 
produces absorption features that evolve out to z > 6. 



5 DISCUSSION 

We have created and analysed predictions for the UV back- 
ground that are intended to broadly span the possibilities in 
star-formation rate and quasar-luminosity density. Our fidu- 
cial model with the HRH07 QSO LF (first entry in Table [TJ 
provides a reasonable match to the level of ionizing fiux in- 
ferred from Lyof forest measurements when an H I escape 
fraction of 0.1 is assumed; this corresponds to a total /esc 
of 1 to 1.5 per cent when combined with the dust absorp- 
tion values predicted by our semi-analytic galaxy formation 
model. The 'low' model for the SFR density, with the lower 
CDM power spectrum normalization of WMAP3, requires a 
larger escape fraction to match ionization rate data, espe- 
cially the higher redshift points z > 5, where /esc, hi ~ 0.5 
is required. Both of these models fail reproduce the nearly 
constant ionization rate between 2 < z < 4.5 seen in some 
fiux decrement analyses. The high-peaked model, which has 
a star-formation rate that increases until redshift 5, does 
produce a somewhat fi a tter io nization rate, as suggested by 
iFaucher-Giguere et al.l (|2008l ). The large amount of high- 
redshift star-formation in this model is not supported by 
estimates of stellar mass buildup, and should be considered 
a somewhat extreme scenario. Another mechanism for pro- 
ducing a fiatter ionization rate history is an evolving es- 
cape fraction that increases with redshift, a possibility that 
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Figure 11. Attenuation factors (e"'') as a function of gamma-ray energy for the indicated source redshifts. Curves are as in Fig. 1101 
and indicate the absorption resulting from our models of star-formation and quasar emissivity. Black solid line: fiducial model with H 
I escape fraction of 0.1. Dashed blue line: low star-formation model, with escape fraction 0.2. Orange dash-dotted line: high-peaked 
star formation rate with escape fraction 0.1. Green long-dashed: fiducial model with SB03 quasar contribution and escape fraxition 0.02. 
Curves for the high-pcaked star formation and high quasar models converge to the fiducial model for z ^ 2. 
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Table 2. Here we summarize a broad set of models, and the qualitative level of agreement of each with Lyo forest data, proximity 
effect measurements, and simulations of total ionizing escape fraction from star-forming galaxies. The numbers of the models shown in 
Figs. [To] [11] and [12] are in bold. 





Model Parameters 








Fits with Data 






Model 


bb H Density 


Quasar Luminosity 


fescHI 


Flux Decrement" 


Proximity h/iiect 


Softness'^ 


J esc 


1.1 


Fiducial 


HRH07 


0.02 


X 


X 


o 




1.2 


Fiducial 


HRH07 


0.10 


o 


X 




%/ 


1.3 


Fiducial 


HRH07 


0.20 


o 


X 


o 




2.1 


Low 


HRH07 


0.10 


X 


X 


o 




2.2 


Low 


HRH07 


0.20 


o 


X 




%/ 


2.3 


Low 


HRH07 


0.50 


o 


X 


o 


o 


3 


Fiducial High-peaked 


HRH07 


0.10 


V 


X 


V 


%/ 


4 


Fiducial 


SB03 Model C 


0.02 


X 


V 


X 





- best agreement 
o - marginal agreement 
X - poor agreement 

° This refers to the relatively flat ionization rate for H I at a level seen in the quasar spectra data of lBolton et al.l ll2005h and 
iFaucher-Giguere et all l|20odb (Fig. [7]|. 

^ The higher levels of ionization rate determined in uncorrected proximity effect measurements (Fig. [7]l. 
" The softness 77 = iV(HeII)/A''(HI), from data compiled in Fig. |1] 

The total escape fraction; recall from Section 2 that this is equivalent to /esc, HI * /l500- In our semi-analytic models, /{^qq ranges 
from about 7 to 14 at the redshifts of interest. 

tant {z 3) galaxies ((Shatalev et al.|[2006l ). but not closer 
sources (ISiana et al.ll2007l ). This would have only a weak ef- 
fect on the opacities we have calculated, as most attenuation 
of gamma rays is due to the non-ionizing UV background, 
which would have a much larger number density than the 
ionizing background even for a high escape fraction. 

Our results suggest that observations of sufficient num- 
bers of high-energy gamma-ray sources out to high red- 
shift could provide a probe of the UV background at these 
epochs that is independent of any other observational test. 
Pair-production with target background photons produces a 
spectral cutoff at energies that are dependent upon redshift 
and assumed cosmological model. With enough detections 
of blazars and/or gamma-ray bursts (GRBs) at different 
confirmed redshifts in the 10 to 100 GeV energy decade, 
it should be possible to detect an evolving cutoff in energy, 
and distinguish between the different background levels pro- 
posed in this work. The exact number of blazars that will 
be detected at GeV energies over the coming years is uncer- 
tain and dependent upon the poorly-understood emission 
processes and number density evolution of these sources. 
However, even conservative estimates indicate that a large 
number of sources will be detectable by the Fermi space- 
craft, which was launched in June 2008 and is currently in 
all-sky survey mode. The EGRET experiment on the Comp- 
ton Gamma-ray Observatory (CGRO) detected more than 
60 high-confidence blazars at energies of > 100 MeV out 
to redshift 2.2 8, mostly of the flat-s pectrum radio quasar 
(FSRQ) type (|Mukheriee et al.lll997D . An extrapolation of 
these results suggests that Fermi will see ~ 1000 blazars 
extending to higher redshift (|Dermerll2007l ). An analysis of 
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Figure 12. The redshifts at which the universe becomes optically 
thick (r > 1) to gamma rays at a given observed energy. Line 
colours and types are as in Fig. 1111 



we do not e xplore here but which has been seen in some 
simulations (|Razoumov fc Sommer- Larsenll2007l . 120061 ) ■ and 
as already mentioned may be suggested by observations 
which have detected Lyman continuum radiation from dis- 
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tw o different realizations of th e blazar luminosity function 
by IChen. Reyes fc Rit j (|2004l ) suggested that Fermi could 
detect thousands of blazars, and would potentially be able 
to measure attenuation at distances as great as 2; = 5. The 
3-month Fermi LAT surve y has already rep orted 106 AGN 
sources at high confidence l|Abdo et al.ll2009l 'l. In addition to 
analyzing blazar spectra in survey and pointed observations, 
Fermi will also act as a finder for new and upcoming ground- 
based experiments such as H.E.S.S.-II and MAGIC-II which 
will be capable of resolving most of the energy ranges of in- 
terest. In survey mode, Fermi will also act as an alert system 
for fiaring sources. 

None of our models predicts significant attenuation at 
10 GeV or below for any redshift. This is true even for our 
extremely quasar-dominated model, where the opacity of a 
10 GeV observed gamma ray is never higher than r ~ 0.2. 
As the ionizing flux in this model is higher than allowed 
by most measurements of the Lya forest, it is unlikely that 
any cosmological model could produce significant gamma- 
ray opacity due to a large contribution of ionizing photons 
to the background. The high-peaked star formation model 
produces the most absorption in the 10-100 GeV energy 
range for z > 3, but despite having a very high UV output 
only has a moderate impact on the calculated optical depths 
relative to the fiducial model. 

5.1 Comparison with other work 

It is useful to compare the absorption predicted by our mod- 
els with the calculations of other authors who have used dif- 
ferent methods, in the cases where their results include our 
energy and redshift regime of interest. In many instances, 
our predicted attenuation is less than what has been previ- 
ously proposed. 

The background model of iFranceschini et all (|2008l l IS 
based upon extrapolated luminosity functions determined 
from a large compilation of multiwavelength data, includ- 
ing deep ACS imaging of distant galaxies, and treats sep- 
arately the evolutionary histories of spiral, elliptical, and 
star-bursting galaxy populations. While the ir EBL agrees 
well w ith our fiducial model at z=0 and z=l (|Gilmore et al.l 
I2OO9I ). their absorption r in the 10-100 GeV energy decade 
is at least a factor of two greater at z =2-4 than any of our 
mode ls. The most recen t mod els of Stecker and collabora- 
tors (|Stecker et al. I I2OO6I . I2OO7I ) are based on a 'backwards 
evolution' model in which galaxies' emission SEDs are de- 
termined by their brightness in one band, taken to be 60 
Hm. The luminosity of the galaxy population at this wave- 
length is assumed to brighten with redshift as a power law in 
{1 + z). One disadvantage of this method is that it attempts 
to describe luminosity evolution over several orders of mag- 
nitude in wavelength from a single broken power law, which 
cannot take into account the complexity of galaxy evolution. 
Gamma-ray opacities in this work are much higher than our 
predictions, with the universe optically thick (r > 1) to 10 
GeV gamma rays above z ~ 3, and for > 25 GeV above 
z — 1. This level of absorption holds very different impli- 
cations for experiments such as Fermi. At high redshifts, 
absorption cutoff spectral features would be visible between 
about 5 and 20 GeV, with no signal from higher energies 
due to optical thickness from the background. The galaxy 
SEDs in these models have no emission above Lyman en- 



ergies, and therefore all attenuation at these very low en- 
ergies is the result of near-threshold interactions with non- 
ionizing UV photons. The redshift-dependent optical and 
U V SEDs used are bas e d on t he population synthesis models 
of lBruzual fc Charlod l|l993h . This model does not include 
UV dust extinction, which as we have found in this work 
can reduce far-UV emissivity by a factor of ~10 at higher 
redshifts. 

The recent observation by Fermi of high-energy emis - 
sion from GRB 080916C at z=4.35 (|Greiner et al.l |2009| ) 
provides a valuable first test of these predictions for GeV 
absorption. The highest energy photon seen by the LAT 
was 13.2 GeV, w ith over 10 photons seen above 1 GeV 
(|Abdo et al.ll2009l ). In all of our models, the gamma-ray op- 
tical depth is much less than 1 for this energy and redshift, 
and similar values are found in the star-formation models of 
Razzaquc ct al. (2009). The models of Stecker and collab- 
orators predict a much higher opacity, r = 3.5 to 4.5, for 
the 13 GeV photon. Including the la error on redshift and 
photon energy in finding maximal and minimal values, the 
corresponding transmission probability could be as high as 
8.2 per cent for their 'baseline' model, or as low as 0.5 per 
cent for the 'fast evolution' model (Stecker, private commu- 
nication). While it is difficult to draw conclusions from a 
single event, more bursts seen with GeV emission equal or 
greater to GRB 080916C could strongly disfavor such a large 
background fiux, and observations of slightly higher energy 
photons (E ~ 30 GeV) from similar redshifts could provide 
a test of our models. 

5.2 Caveats and future work 

We expect our approach to be reasonably accurate at pre- 
dicting the ionizing and non-ionizing background fields out 
to redshift ~ 6, where H I Gunn- Peterson troughs appear 
in observed quasar spectra (|Fan et al.ll2006l l. At higher red- 
shifts, during the epoch of reionization, the concept of a uni- 
form background for ionizing photons is no longer valid, as 
photons above the Lyman limit are confined to the vicinity 
of their sources. In our Lyman absorption model (Equation 
[3|, the sudden increase seen in H I opacity at redshift six is 
not represented, and our model would therefore be expected 
to overproduce the ionizing background above this redshift. 
A similar limitation exists in our treatment of He opacities 
above the redshift of He reionization z ~ 3. However, these 
factors alone are unlikely to have a significant effect on cal- 
culated opacities. 

We have made the assumption of a universal stellar 
IMF in this work, and have not included a separate pop- 
ulation of metal-free (population-Ill ) stars or other earl y 
source types such as miniquasars (|Madau et al.l l2004al ). 
These types of unobserved sources could have very different 
spectra than standard stellar populations, and could pro- 
duce large contributions to the ionizing and non-ionizing UV 
backgrounds. Because of redshift effects, gamma rays with 
low (< 10 GeV) observed energies for very high-redshift 
sources could have significant interactions with the freely- 
propagating non-ionizing background. It is therefore possi- 
ble that opacities at reionization redshifts could be much 
higher than we propose here due to unseen UV production 
mechanisms. While models for gamma-ray blazars do not 
typically predict sources at these high redshifts, GRBs are 
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known to exist above redshift six (jCreiner et al.ll2"009al ) . and 
long-duration GRBs could potentially be seen as far out as 
star-formation occurs. The EGRET experiment, operating 
from 30 MeV up to ^30 GeV, was able to view a small num- 
ber of photons from GRBs, and the detection of high-energy 
emission from GRB 080916C by the Fermi LAT demon- 
strates the ability of this instrument to detect GeV photons 
from these events. Though predictions are highly uncertain, 
it is possible that GRBs could produce significant numbers 
of photons well ab ove 10 GeV through inverse-Compton or 
hadro nic processes (|Le fc Dermerll2008l : lAndo. Nakar fc Saril 
I2OO8I ). Calculations of the background flux from some pos- 
sible reionization scenarios and source types at z > 6 may 
therefore be a worthwhile undertaking. 

It has been suggested by a number of authors that the 
discrepancy between observed stellar mass density and in- 
stantaneous star-formation rate density (see Section 12. 2[) 
could be explained by an IMF that evolves with redshift 
or is more top-heavy in r apidly star-forming galaxies (jPavel 
I2OO8I : iFardal et al.ll2007l : iBaugh et aLlbOOSl ). Alternatively, 
an IMF with shallower high-end slope has been suggested 
as a source of early reionization (Chary 2003) • Altering the 
high-mass end of the IMF will change the spectrum pro- 
duced by galaxies and also the attenuation by dust, although 
as probes of star-formation generally involve the same high- 
mass stars that produce the UV background, there is some 
degree of degeneracy between these two quantities when the 
assumed IMF is changed. This is also an issue that warrants 
further study. 
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